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. ABSTRACT 

' Using the IRAM Plateau-de-Bure Interferometer (PdBI) we have searched for the 

, upper fine structure fine of neutral carbon (Ci(^P2 ^Pi), t'rcst = 809 GHz) and 

00 ■ ^2co(J=7^6) (i/rest=806GHz) towards the submiUimetre galaxies (SMGs) GN20 

O ; (SMMJ123711.9+622212, z = 4.055) and GN20.2 (SMM J123708.8+622202, z = 

0^ . 4.051). The far-infrared (FIR) continuum is detected at 8f7 significance in GN20, with a 

■ flux density of Si.smm = 1.9±0.2mJy, while no continuum is detected in GN20.2. Both 
^ I sources are statistically undetected in both Ci(^P2^ ^Pi) and ^^CO(J=7^6) lines; we 

■ derive line luminosity limits for both Ci and CO of L' <, 2 x 10^" Kkms^^ pc^. As- 
' suming carbon excitation temperatures of rex = 30K (the galaxies' measured dust 

^ ! temperatures), we infer Ci mass limits of Mci < 5.4 x lO*^ M0(GN2O) and Mqi < 

■ - - ■ 6.8 X 1O6M0(GN2O.2). The derived Ci abundance limits are < 1.8 x IQ-^ for GN20 

and < 3.8 x 10~^ for GN20.2 implying that the systems have Milky Way level neutral 
carbon enrichment (X[C i]/X[H2]) or lower, similar to high-redshift carbon-detected 
systems (at 5 x 10~^) but about 50 times less than the neutral carbon enrichment of 
local starburst galaxies. Observations of GN20 and GN20.2 in high-resolution MER- 
LIN-I-VLA radio maps of GOODS-N are used to further constrain the sizes and loca- 
tions of active regions. We conclude that the physical gas properties of young rapidly 
evolving systems like GN20 and GN20.2 are likely significantly different than star- 
burst /ULIRG environments in the local Universe yet similar to z 2 SMGs. Unless 
gravitationally amplified examples can be found, observations of galaxies like GN20 
will require the order of magnitude increase in sensitivity of the Atacama Large Mil- 
limetre Array (ALMA) to constrain their Ci and high- J CO content, despite the fact 
that they are the brightest systems at z ~4. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: infrared — 
galaxies: starbursts — galaxies: individual (GN20/GN20.2) 
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1 INTRODUCTION 

Examining the molecular gas content of distant galaxies 
is fundamental to our understanding of galaxy formation 
and evolution theories. Recent observations have shown that 
many distant objects contain gi ant molecular gas reser- 
voirs (> lO'^" Mq; for a review see lSolomon fc Vanden Bout! 
I2OO5I ). These gas repositories are thought to provide 
the fuel needed to power the extreme starbursts ob- 





e.g. iFraver et al.l 


19991: iNeri et al. 


200a 


:lGreve et al.ll2005l:lTacconi et al.ll2006 


: ChaDman et al. 



2008 



In addition, detection of AGN in molecular gas sug- 
gests a link between the most massive starbursts, the growth 
of mas sive black holes an d the onset of strong nuclear ac- 
tivity l|Coppin et al.ll2008l ). Molecular gas has already been 
detected in about 60 high redshift sources, fro m star forming 
galaxies to quasars, at redshifts 2 < z < 6 .4 ("Walter e t al. 
2003 ; Bcrt oldi et al. 2003; Grov e et all [200 5: Tacconi e t al. 
20061 . 120081 ). 



Carbon monoxide has bright rotational transitions and 
is therefore the most commonly observed tracer of molec- 
ular gas; however, CO lines are usually optically thick and 
difficult to model. In contrast, observations of optically thin 
neutral carbon (C i) can be used to derive gas physical prop- 
erties without requiring detailed radiative transfer models. 
This is because carbon has a 3P fine-structure forming a 
simple, easily-analyzed three-level system. The gas excita- 
tion temperature, neutral carbon column density and mass 
may be derived independent of any other information pro- 
vided there are detections of both carbon lines, C ij^Pi — > ^Pq) 
(492 GHz) and Ci(^P2-* ^Pi) (809 GHz) (e.g. IStutzki et all 
ll997l : IWeiss et al.ll2003l . l2005l ). Furthermore, observations of 
either Ci line could potentially be used as a probe of in- 
ternal gas distribution which is independent from the more 
luminous optically thick CO observations. 

While the Earth's atmosphere has very low transmis- 
sion at the rest frequency of Ci(''Pi ■^Po) and is virtually 
opaque at Ci(^P2 ^Pi), some studies of neutral carbon 
have been performed in nearby sources where the lines are 
bright enough to be observed. These include the Galactic 
Centre, molecular clo uds in the galactic disk, M82 and other 



nearby galaxies (|e.g. White et al.l|l994j : jStutzki et al.lll997l: 



^, this sug- 
0) transitions 



Gerin & Phillips"2000;'Oiha ot al.'200l': Is rael fc Baasll2002 : 
Schneider et al. 2003; Kramer ot al. 200j). These studies 
show that C I and CO emission trace each other well, in- 
dependent of the type of heating environment. Since they 
share similar critical densities, Ud ~ 10^ cm 
gests that the Ci{^Pi^ ^Po) and ^^C0(J=1 
arise from the sam e gas volume and have similar excita- 
tion temperatures (jlkeda et al.ll200^ ). In addition, several 
studies have found excellent agreement bet ween C i and CO 
derived H2 masses in local ULIRG s (e.g. |g erm fc Phillips! 
ll998l : |Papadopoulos fc Grevell2004l ). 

Despite improved atmospheric observing conditions at 
the redshifted (z > 2) C i frequencies, observations of dis- 
tant, faint sources are particularly difficult; neutral carbon 
has only been confirmed prev iously in five other high red- 
shift sources, th ree at z ^2.5 JWeiss et al.ll2003l . l2005h . one 
z = 3.91 QSO (IWagg et"all I2OO6I ) and one z = 4.12 QSO 
l|Petv et al. 12004 ). While it is challenging to detect, Ci is a 
cooling line and therefore is important in understanding the 



composition of a galaxy's dense interstellar medium (ISM). 
In high redshift sources, detection of C i indicates that the 
ISM has condensed and become significantly enriched while 
the Universe was still very young. 

Here we report on the search for the upper fine struc- 
ture line of neutral carbon, Ci('^P2 — > ^Pi), as well as the 
^■^C0( J=7— >6) line, towards two of the highest redshift and 
most luminous submillimetre galaxies, GN20 {z = 4.055, 
SMMJ123711. 9+622212, identified in Pope etal., 2006 
whose redshift was accurately measured in iDaddi et al.l 
|2009| . hereafter referred to as D09) and GN20.2 {z = 4.051, 
SMM J123708.8+622202, also detect ed in ^^CO( 7=4-^3) i n 
D09; originally from the catalogue of IChapman et al.ll200ll ). 
GN20 (5850=20.3 mJy) and GN20.2 (5850=9.9mJy) are two 
of th e brightest submillimetre galaxies in GOODS-N (see 
also |Popdl2007h . Both of their rest-frame ultraviolet spec- 
tra lack any emission features (D09) and a re similar to the 
UV s pectra of many other z ~ 2.5 SMGs (jChapman et al.l 
120051) . The absence of AGN characteristics in optical spec- 
tra s uggests that the g as is mainly heated by star forma- 
tion. IPope et al.l (|2006l ) have also shown that GN20 has a 
very similar SED and mid-IR spectral properties to other 
SMGs, as are its X-r ay luminosity and photon index (cf. 
[Alexander et al.ll2005l ). 

In addition to our interferometric millimetre observa- 
tions from the Plateau de Bure Interferometer, we present 
high- resolution radio maps of both galaxies at 1.4 GHz from 
MERLIN-I-VLA to further characterize these bright SMGs. 
We organize the paper as follows: observations are described 
in Section (2] results are discussed in Section O and our 
discussion of the implications on high-redshift SMC enrich- 
ment are given in Section 3] Throughout, we use a A CDM 
cosmology with Hp = 71kms~^ Mpc"'^ and ilm = 0.27 
(|Hinshaw et al.ll2009l) . 



2 OBSERVATIONS 

2.1 Molecular Line Observations 

Observations were carried out with the IRAM Plateau- 
de-Bure interferometer through August 2008, with a 5 
dish D-configuration (i.e. compact). We used the 2-mm 
receivers tuned to 159.873 GHz, midway between the ex- 
pected redshifted frequencies of the Ci(^P2 — > ^Pi) {vieat = 
809.342 GHz) and ^'^CO{J=7-^6) (f^cst = 806.651 GHz) 
transitions for GN20 at 2 = 4.055 and GN20.2 at z ^ 4.051. 
The pointing centre was closer to GN20 than to GN20.2; 
since they are both off phase center, their fluxes require pri- 
mary beam attenuation correction factors of 1.11 (GN20) 
and 2.62 (GN20.2). This observation was made possible be- 
cause of the redshift measurement from the ^■^CO( J=4— >3) 
D09 observation. The synthesized beam size at 160 GHz 
is ~3". Calibration was obtained every 12min using the 
standard hot /cold-load absorber measurements. The source 
3C454.3 was used for absolute flux calibration. The antenna 
gain was found to be consistent with the standard value of 
29 Jy K^^ at 160 GHz. We estimate the flux density scale to 
be accurate to about ±15%. 

Data were recorded using both polarizations overlap- 
ping, covering a 900 MHz bandwidth. The total on-source 
integration time was Shrs. The data were processed using 
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Table 1. Observed and Derived Properties of GN20 and GN20.2 





01N2U 




f^-^MERLIN+VLA (J2000) 


12:37:11.96 


12:37:08.80 


DeCMERLIN+VLA (J2000) 


+62:22:12.4 


+62:22:01.9 


RA160 (J2000) 


12:37:11.89 




Decieo (J2000) 


+62:22:12.4 




(MJy) 


73.8+12.8 


170.0+12.8 


Ssmm" (mjy) 


0.33 


<0.2 


•S'l.2mm'' (mJy) 


9.3 


<2.7 


Si.imm" (mJy) 


11.5 




5850" (mJy) 


20.3 


9.9 


Ri/2 MERLIN+VLA (") 


0.38+0.15 


0.30 


13 


1.4 (derived) 


1.5 (fixed) 


Lfir (L©) 


l.OxlO" 


5.0x10^2 


Tdust (K) 


30+4 


30+12 


^CO[4-3l'' 


4.055+0.001 


4.051+0.003 


-fcOK-S]' (Jykms-1) 


1.5+0.2 


0.9+0.3 


^CO[4-3]' (Kkms-lpc2) 


6.2x101^^ 


3.7x10^° 


Mh,"''' (Mq) 


5.0x101" 


3.0x10^° 


Sieo Continuum (mJy) 


1.9+0.2 


0.5+0.3 


-''co[7-6] =-fci (Jykms-l)'* 


<1.2 


<1.9 


-^00(7-6] =-^01 (Kkms-lpc2) 


< 1.6 X 10"' 


< 2.5 x 10^° 


(from CO[7-6]) (Mq) 
MciIto, (Mq) 


< 1.3 X lO^" 


< 2.0 x 10" 


< 5.4 X 10^ 


< 6.8 X 10^ 


X[CI]/X[H2] 


< 1.8 X 10"^ 


< 3.8 X 10"* 



Table Notes. We fit modified blackbodies to the FIR continuum 
fiux densities (see section 13. ip to derive dust temperature and 
^FIR- is reliably constrained to 1.4 for GN20 and fixed to 1.5 
for GN20.2. 

° Obse rved continuum fiux densi ties from the liter ature at 
850um dPope et a.l.ll2006h . 1.1mm jPerera et al.ll2008l ). 1.2mm 
jGreve et al.ll2008l) . and 3mm (D09). Upper limits for GN20.2 
are 2a. The radio fiux (5i.4) is measured at 1.4 GHz from VLA. 
' Observed properties from D09, given here for comparison. 

Derived from the measured D09 CO [4-3] line strength. 
D09 assumes a constant brightness temperature between 
12CO(J=1^0) and 1^00(7=4^3) to derive Mh2, but 

= 8/16 (cf. 

"CO[7-6]'-"CO[l-0] ^ ^/"^^ 

IWeiss et al.l 12007^ . We also assume a con version factor X = 
M^/L^Q=0.8 M0 (K km s-i pc^)-! (cf. lDownes fc SolomonI 
1 199^ . 

Line intensity limits are 2a. In this data, both Ico a-nd Id 
limits are equal since the noise measurement in each frequency 
range is equal, as are the assumed line widths (700 kms^^). 



here we assume that ^co[4— 3]/^co[l 
IWeiss et al] l2007t) . and L' JL' 



the GILDAS packages CLIC and MAPPING and analyzed 
with our own IDL-based routines. The extracted data cube 
(2 sky coordinate axes and 1 spectral axis) has an RMS noise 
of 0.6 mJy. For clarity of presentation, we have re-gridded 
the data to a velocity resolution of ~37 kms~i (20 MHz). 
No obvious {>5a) emission lines are seen in the data; how- 
ever, the next section statistically tests for the presence of 
FIR continuum and CO and C i emission lines. 



2.2 MERLIN+VLA Radio Imaging 

High-resolution observations from the Multi- Element Ra- 
dio L inked Interferometer Network (MERLIN; IXhomassonl 
Il98d)_ were obtained for these sources as described 
in IMuxIow et al.l (|2005l ). with an RMS noise of ~ 
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Figure 1. MERLIN+VLA radio contour overlay on H ST /ACS 
BVi colour images of GN20 and GN20.2. The field sizes are 
3" X 3" and the MERLIN beam size is 0.3" (outer contour of 
corner inset is the FWHM of the beam). The levels of the ra- 
dio contours are drawn at 3, 4, 5, 7, and 10 cr; the RMS noise in 
the MERLIN+VLA images is 2.9 /ijy and 4.6 ^tjy for GN20 and 
GN20.2 respecti vely. The centre o f the SM A continuum emis- 
sion analyzed by llono et"aL I ||2006| ') and [Youneer et al. j ||2008|) is 
marked by the green cross on the GN20 map. The emission cen- 
troid for GN20 in the mid-IR {Spitzer IRAC/MIPS) is consistent 
with the radio and FIR continuum position. The FIR-optical off- 
set for GN20 is significant at the 3 a level. 



5^tJybe am~i. While GN20 o r GN20.2 were not imaged di- 
rectly in lMuxlow et all Hooi), the source is within the sen- 
sitive region of the data (Lovell Telescope primary beam) . A 
combined MERLIN+VLA map was then constructed, with 
a sensitivity of 3-4 ^Jy beam~i and high positional accuracy 
(tens of mas). In Figure [l] we show MERLI N radio contours 
on top of the HST BVi tricolor images (|Giavalisco et al.] 
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Figure 2. The SED fits to the FIR measured photometry of 
GN20 (black) and GN20.2 (gray). The inferred continuum flux 
densities at 1.88 mm are shown as large diamonds: 1.8±0.4mjy 
for GN20 and 0.8±0.4mjy for GN20.2. Their full SEDs from 
10 /im near-IR to 10cm radio (rest wavelength) are shown in the 
inset; the area outlined in the dotted box is the zoomed-in region 
shown in the larger plot. 

l2004l fl. The field size is 3" x 3" and the restoring beam 
size is 0.3". We measure the center of MERLIN+VLA ra- 
dio emission as 12:37:11.96, +62:22:12.4 for GN20 and 
12:37:08.80, +62:22:01.9 for GN20.2. 



3 RESULTS 

3.1 The Continuum Contribution at 1.88mm 

Interpreting the detectability and significance of emission 
lines in the millimetre requires an estimation of the expected 
continuum contribution at the wavelength of observations 
(160 GHz, which is a wavelength of ~1.88mm). Indepen- 
dent of our PdBI observations, we estimate the continuum 
flux at 1.88 mm by fitting modified-blackbody SEDs to the 
galaxies' FIR flux measurements. The continuum flux den- 
sities for GN20 and GN20.2 are given in Table 1. We fixed 
the emissivities of the galaxies to (3 = 1.4 for GN20 (see 
|Popa.2007 ) and (3 = 1.5 for GN20.2 is common practice for 
galaxies in which /3 is not well constrained, although we note 
that letting P vary has a minimal effect, <10%, on the con- 
tinuum flux). With these SED constraints we measure the 
dust temperatures and FIR luminosities that are given in 
Table 1. The SED fits to both galaxies may be seen in Fig.[51 
From these flts we infer that the 1.88 mm flux generated by 
FIR continuum should be 1.8+0.4 mjy and 0.8+0.4 mjy for 
GN20 and GN20.2, respectively. 

Having inferred the expected continuum level we must 
now measure the continuum from the 160 GHz PdBI obser- 
vations. A 1-D spectrum for each object is flrst extracted 
at the peak integrated flux position, shown in Fig. |3l If we 
ignore the possible presence of emission lines, the average 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, and obtained from the Hubble Legacy Archive, 
which is a collaboration between the Space Telescope Science In- 
stitute (STScI/NASA), the Space Telescope European Coordinat- 
ing Facility (STECF/ESA) and the Canadian Astronomy Data 
Centre (CADC/NRC/CSA). 



flux over the whole bandwidth is 1.9+0.2 mJy for GN20 and 
0.5+0.3 mJy for GN20.2. The latter is consistent with the ex- 
pected continuum flux, 0.8+0.4 mJy, for GN20.2, therefore 
we conclude that neither Ci(3P2 ^ ^Pi) nor ^^CO(J=7^6) 
is detected in GN20.2. The GN20 average flux is consistent 
with the expected continuum level, 1.8 mJy, but since there 
is a large uncertainty in the expected flux at 1.88 mm, it 
needs to be further examined to see if the flux excess is due 
to the partial detection of C i and CO lines. 



3.2 Placing limits on Ci(3P2^ ^Pi) and 
i2co(J=7^6) 

Figure |3] shows extracted 1-D spectra alongside integrated 
channel maps for GN20 and GN20.2. The maps are centred 
on each galaxy's VLA radio position, and the 1-D spectra 
are extracted at the point of peak integrated flux for GN20 
and at the VLA centroid for GN20.2. Our observations have 
a ~0.6mJy RMS noise near GN20 and 1.6 mJy RMS noise 
near GN20.2across the entire 900 MHz bandwidth. 

To measure the continuum without allowing contami- 
nation from possible line emission, we must extract the por- 
tion of the GN20 spectrum where line emission is antici- 
pated. We use the measured D09 redshift and line width of 
^^CO(J=4^3) (2 = 4.055 and Av ~ 700 kms"^) as an a 
priori condition and we apply them to both ^^CO(J=7— >6) 
and Ci('^P2 — ► ^Pi) lines. Unfortunately, the narrow band- 
width limits the "non-line" spectral portion to a small 
number of channels, speciflcally those with frequencies 
159.75 < z^obs < 159.95 GHz or i/obs > 160.25 GHz. We mea- 
sure the continuum flux density in this spectral region to 
be 5*160 (GN20) = 1.4+0.3 mJy. This suggests that a flux ex- 
cess does exist around the expected lines, from 159.95 < Vahs 
< 160.25 GHz (the C i line region) and at i^obs < 159.75 GHz 
(the CO line region). To test the significance of these ex- 
cesses, we measure the mean flux density in each region and 
compare it to the "non-line" flux. The average flux density 
is 1.6+0.3 mJy in the Ci line region and 2.5+0.4 mJy in the 
CO region. Both excesses are of low signiflcance (<2cr); thus, 
we conclude that neither ^^CO(J=7^6) nor Ci(3p2^ ^Pi) 
is detected in GN20. 

With no signiflcant detection of CO or C i lines, we at- 
tribute all flux across the band to the FIR continuum, there- 
fore the measured continuum flux densities are 1.9+0.2 mJy 
for GN20 and 0.5+0.3 mJy for GN20.2. The latter is con- 
sistent with no detected continuum. The former continuum 
measurement for GN20 is signiflcant at the 8ct level. Al- 
though some flux from the emission lines might contribute 
to the flux density measurement of the continuum, we have 
not accounted for this additional uncertainty in our mea- 
surement. 

We derive 2(j upper limits to C i and CO line intensi- 
ties assuming the redshift and line width given by D09. For 
GN20, we measure Ici ~ Ico < 1.2 Jykms^^ (both limits 
are equal since the assumed line widths and the noise mea- 
surements in each frequency range are equal). For GN20.2, 
we measure Ici ~ < 1-9 Jy km s~^ . The corresponding limits 
to line luminosity are given in Table [T] 
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Figure 3. The 160 GHz maps and spectra for GN20 (left) and GN20.2 (right). The 20" X 20"maps are integrated over the whole 
bandwidth; 4, 5, 6, and 7 cr contours are over-plotted (the detection of GN20 is Scr). The ~3" beam shape is shown as a lower right 
insert in each map. The 20 MHz spectra are extracted at the peak flux position of GN20 and at the VLA centroid position for GN20.2. 
The average flux, or continuum estimation in the absence of line emission, is shown as a horizontal red line bounded by Icr uncertainties 
(red dotted lines). The expected frequency centres of the Cl and CO emission lines (based on the D09 redshift) are marked by dashed 
gray vertical lines. 



3.3 CO Excitation and Brightness Temperature 

For interpreting CO line intensities, the assumption of con- 
stant brightness temperature has been widely taken a priori 
in the past due to lack of relevant data. However, many re- 
cent observations of multiple CO transitions in high red- 
shift sources have shown it to be an unwarranted assump- 
tion (|Dannerbauer et al.ll2009l ). Although detailed radiative 
transfer models demonstrate that the ^ CO(J=1^0) and 
^^CO( J=3— i-2) transitions are near thermal equilibrium (i.e. 
L'rnn-M ~ _pr^n(^ -9)) for the three z ~ 2.5 sources studied 
jWeiss et al.l |200^), there are several AGN dominated galax- 
ies which do not have constant brightness temperature for 
high-J transitions (e.g. FSC 10214-^4724, APM08279-H5255, 
and the Cloverleaf Quasar). Here we adopt the bright- 
ness tempe r ature conversions for SMGs observed by 
IWeiss et al.l (|2007l ). see their Figure 3; explicitly, we as- 
sume flux density ratios of Sco[4-3]/Sco[i-o] = 10.0±0.8 and 
Sco[7-6]/Sco[i-o] =8±1. This implies Lco[7-6]/-t'co[i-o] = 
8/49 and -£'co[7-6]/-'^co[4-3] = 64/245. Working back- 
wards from the D09 ^^CO(J=4^3) detection, this im- 
plies that ^■^CO( J=7— >6) in GN20 would have luminosity 
~ 1.6 X 10^° K km s-i pc^ We measure a line luminosity limit 
of < 1.6 X 10^" Kkms~^ pc^. While there are large uncer- 
tainties in our brightness temperature conversion assump- 



tion, this measured result for GN20 agrees with the esti- 
mated population of the higher-J transitions according to 
other high- 2: SMGs. Our data for GN20.2 are less constrain- 
ing due to the greater noise; our limit on line luminos- 
ity is consistent with ico[7-6] ~ ^-2 ^ 10^" Kkms~^ pc^, 
which is the expected luminosity extrapolated from the 
^^CO(J=4^3) detection. 



3.4 Radio Source Sizes 

We use high-resolution MERLIN-I-VLA radio maps to con- 
strain the size of the 1.4 GHz emission, to better assess the 
contributions from AGN and starbursts. The measurements 
of angular size were made using AIPS software JMFIT, de- 
convolving the restoring beam before attempting to derive 
the fitted Gaussian size. Since the signal to noise is quite 
low in GN20, it is impossible to detect component exten- 
sions smaller than the width of the restoring beam. Figure [1] 
shows the MERLIN-I-VLA 1.4GHz radio maps as contours 
on top of HST ACS i-band imaging. We measure the angu- 
lar FWHM of the GN20 emission region to be 0.38±0.15", 
but it is only significant above the 0.3" beam size at the 
2.4(T level, while the 0.5" radio-optical offset is significant at 
the 3(7 level. 
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Ilono et al.l (|2006l ') showed that the 850 //m Smithsonian 
Millimeter Array (SMA) position of GN20 was significantly 
offset from the optic al emission (0.8" offset also at the So- 
level) . I Younger et"all (fioOSi ') obtained higher resolution con- 
tinuum imaging of GN20 with the SMA and constrained 
the size to 0.60 ± 0.13" claiming it to be partially resolved. 
This is consistent with our MERLlN+VLA-measured size 
of 0.38±0.15" due to the large uncertainties on both mea- 
surements. Because GN20 is not exceptionally bright in the 
radio (74 /^Jy), we cannot determine if its asymmetric, ex- 
tended MERLIN-I-VLA morphology is significant; therefore, 
we cannot rule out that GN20 might be a radio point source 
dominated by AGN activity. 

The radio source size of GN20.2 is measured to be en- 
tirely unresolved at < 0.30". Radio emission with 51.4 > 
100 ^Jy at 2 = 4.051 constrained by R\/2 < 0.30" would 
have to emit at a super-Eddington rate if it were completely 
dominated by star formation (a calculation ba sed on the 
theoretical maximum st ar formation density, cf. lElmegre"enl 
1 19991 : ICasev et allbOOgl ). For this reason, GN20.2's radio 
emission is likely partially generated, if not dominated by 
an AGN at the galajcy's core. 



4 DISCUSSION 

4.1 Neutral Carbon and H2 Mass 



Using the formulae given in IWeiss et al.1 l|2003h for the upper 
fine structure line of neutral carbon Ci('^P2 ^Pi), we can 
derive the mass of neutral carbon via 



Mci = C mci 



hc^Aio 



CI(3P2^3pj), 



(1) 



where ^(Tex) = l-|-3e"^i''^'=- -l-Se'^^/^"- is the neutral car- 
bon partition functior0. In the equation above, Ti = 23.6 K 
and r2=62.5K represent the energies above the ground 
state, L'qi(^3p^_^ 3p_^^ is given in K kms~^ pc^, mci is the mass 
of a single carbon atom, and C is the conversion between cm^ 
and pc^, 9.5x10^® cm^^/pc^. This assumes both carbon lines 
are in local thermodynamical equilibrium (LTE) and that 
C I emission is optically thin. Inserting numerical values in 
Equation [T] yields 



Mci = 4.566X 10"" Q(Ta^) - e 



Lci(3p,^3p,)[Mq].(2) 



Determining Mci definitively requires measurement of the 
excitation temperature, which can only be done by detecting 
both C i(3Pi ^Po) and C i{^P2 ^Pi), the two C i fine struc- 
ture lines. If Tex cannot be measured directly, it is often as- 
sumed to be equal to the characteristic dust temperature (al- 
though in the case of the Cloverleaf Quasar, Weiss et al. 2003 
showed that Tox^SOK, significantly less than the dust tem- 
perature, 50 K). As the excitation temperature is unknown 
for both GN20 and GN20.2, we assume that Tex = Tdust, 
which is appropriate, since here we are only able to derive 
upper limits on carbon mass for both systems. If the exci- 
tation temperature is less than the dust temperatures (like 



^ We note that the negative sign in the ex ponent of 1/5 e 
was accidentally o mitted from the text of iWeiss et al.l ||2003|) and 



the Cloverleaf), the upper limits on carbon mass would de- 
crease. Wit hout account ing for the possible magnification 
proposed bv iPopel l|2007l ). we derive upper limits on carbon 
mass of 5.4 x 10*^ Mq (GN20) and 6.8 x 10'' M© (GN20.2). 

Table [1] also lists the limits on H2 mass derived 
from the ^^CO(J=7— >6) line luminosity li mits, assum- 
ing b rightness temperature conversions from IWeiss et al.l 
({ 2007' ). We also assume the standard ULIRG conver- 
sion factor X = Mu^/L r^c, = 0.8 M© (Kkms"^ pc^)"^ 
(iDownes fc Solomon|[l998l ). 



4.2 C I Abundance in the Molecular Gas 
Reservoir 

Neutral carbon abundance, which is independent of 
magnification and adopted cosmology, is estimated us- 
ing the ratio of masses between C i and H2 through 
X[CI]/X[H2] = Afci/(6MH2). We calculate acarbon abun- 
dance upper limit using the value derived by D09. 
We find that X[CI]/X[H2] ;S 1.8 x 10"^ for GN20 and 
X[CI]/X[H2] ;S 3.8 X lO"'' for GN20.2. Both limits are con- 
sistent with the Ci abunda nce in the Milky Wa y Galaxy, 
A:[CI]/X[H2] = 2.2 X 10~^ (jFrerking et al.lil989l ). but are 
much lower than the measured carbon abundances of local 
starburst galaxies ~ 1.5 x 10~'^ (~50 t imes the abundanc e 



IWeiss et al.l ^m^. 



of the MW, e.g. ISchilke et al.l Il993l: IWhite et al.l Il994l l. 
Ilsrael fc Ba"a3 l|200ll) and llsrael fc Baad (|2003l ') find similar 
carbon abundances in local ULIRGs although they account 
for total carbon mass instead of C i alone. In comparison, 
the handful of carbon abundance measurements that have 
been made at high redshift averag e to ~ 5 ± 3 x 10~® (e.g. 
IWeiss et af]l2005l : IPetv et al.ll2004h . This is consistent with 
both the Milky Way's abundance and the two limits we have 
derived for GN20 and GN20.2. We note that an abundance 
near (or exceeding) that of the Milky Way implies that the 
cold molecular gas in GN20 or GN20.2 would already be sig- 
nificantly enriched 1.6 Gyrs after the Big Bang (a possibility 
that is not excluded by our limits). More observations of C i 
and CO in high redshift galaxies and QSOs are needed to 
investigate possible C i abundance variations with redshift 
or variations between differently heated environments (star 
formation dominant versus AGN). 



5 CONCLUSIONS 

GN20 has become one of the best studied high redshift star- 
burst systems in the GOODS-N field, due to its high ap- 
parent infrared luminosity. Its diverse multiwavelength cov- 
erage has revealed very low apparent levels of AGN emis- 
sion relative to its high star formation. Observations of 
GN20 in molecular emission lines have allowed a charac- 
terization of the system's gas properties, from the strong 
detection of ^^CO( J=4^3) of D09 to the non-detections of 
Ci(3p2 ^ ^Pi) and '2CO(J=7-^6) in this paper. We have 
also characterized the molecular gas of GN20.2, a nearby 
SMG at = 4.051. 

Millimetre and Sub-millimetre continuum observa- 
tions suggest a continuum fiux density at 160 GHz of 
5'i60GHz = 1.8mJy for GN20 and 0.8 mjy for GN20.2. We 
measure 1.9±0.2 mJy continuum for GN20 and 0.5±0.3 mJy 
continuum emission from GN20.2. The former is significant 
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at the 8a level while the latter is not statistically different 
from zero. 

GN20 and GN20.2 are undetected in C i{^P2 ^Pi) and 
^^CO(J=7^6) using observations from the IRAM-PdBI in 
the D-configuration at 160 GHz. Line intensity upper limits 
are given as 7co[7~6](GN20)= 7ci(GN20) < 1.2Jykms-\ 
and /co[7-6](GN20.2)= /ci(GN20.2) < 1.9Jykms"\ 

High-resolution radio imaging from MERLIN+VLA 
show that the radio emission of GN20 is extended over 
0.38±0.15" (2.7 kpc). This size i s consisten t with the 
FIR continuum size measured by llono et al.1 l|2006l) and 
lYounger et al.l (|2008l '). In contrast, GN20.2 is unresolved at 
<0.30" (2.0 kpc), suggesting that AGN emission may dom- 
inate its radio emission. 

The conditions for metal enrichment of the ISM of 
GN20 and GN20.2 do not appear to differ greatly with other 
high redshift sources or the Milky Way. However, local star- 
burst galaxies have ~50 times the carbon abundance of the 
high-z systems and the upper limits for GN20 and GN20.2. 

Future observations from the Atacama Large Millimeter 
Array (ALMA) of cooling gas emission lines like C i{^Pi 
^Po) and Ci(^P2 — » ^Pi) in high redshift sources will signif- 
icantly advance the understanding of the relationship be- 
tween early star formation and metal enrichment. Although 
difficult (in the absence of a boost from gravitational lens- 
ing) to observe with current instruments as evidenced by 
the lack of detection in GN20, one of the brightest known 
high redshift SMGs, C i is an excellent tracer of cold molec- 
ular gas in galaxies. It may even provide useful constraints 
on gas content when undetected relative to CO, and may 
be used as a powerful diagnostic of galaxy evolution out to 
the highest observable redshifts, especially in the light of the 
next generation of mm and submm telescopes. 
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